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e Octave is the "open-source Matlab" 
e Octave is a great gnuplot wrapper 
e www.octave.org 


e www.mathworks.com 


Octave and Matlab are both, high-level languages and mathematical 
programming environments for: 


e Visualization 
e Programming, algorithm development, prototyping 


e Scientific computing: linear algebra, optimization, control, statistics, signal 
and image processing, etc. 


Beware: Octave/Matlab programs can be slow 
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Comparison Matlab vs. Octave 
e Matlab is more flexible/advanced/powerful/costly 

e Has an IDE and a proper GUI, windows for help, variables, command history, etc. 
e Octave is for free (GPL license) 


e There are minor differences in syntax 


This tutorial 
e Applies to Matlab AND Octave unless stated otherwise! 


e |s valid for the 2009 versions 
e Octave 3.2.3 
e Matlab 7.6 


or higher 
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e To start Octave type the shell command octave or whatever your OS 
needs. 
You should see the prompt: 


octave:l» 


e Matlab will start its own window-based development environment 


e If you get into trouble, you can interrupt Octave by typing Ctr1-C 


e To exit Octave, type quit or exit 


O i Arras 


Start, Quit, Getting Help ee es 


* To get help, type help or doc 


。 To get help on a specific command (= built-in function), type 
help command 


e Examples: help size, help plot, help figure, 
help inv, … 


* To get help on the help system, type help help 


e In Octave: type q to exit help mode (like man pages) 


Social Robotics Lab 
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In the help text of Matlab functions, function names and variables are in 
capital letters 


Example: help round returns 


ROUND Round towards nearest integer. 


ROUND(X) rounds the elements of X to the nearest 
integers. 


See also floor, ceil, fix. 


[asa] 


Don't get confused! The naming convention specifies lowercase letters 
for built-in commands. It is just a way to highlight text 


Octave texts are mixed, in lower- and uppercase 


Update: this finally changed in new Matlab versions 
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In Octave/Matlab almost everything is a matrix! 
e Matrices 


(“Matlab” comes from "Matrix Laboratory”) 


Main matrix classes 
e Strings: matrices of characters 
e Structures: matrices of named fields for data of varying types and sizes 


e Logical: matrices of boolean 0/1-values 


Not treated in this tutorial 
e Cells (like structures) 


e Function handles (pointer to functions) 
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What about... 


e Vectors or arrays? 
— A matrix with one column or row 


e Scalars? 
— A matrix of dimension 1x1 


e Integers? 
— A double (you never have to worry) 


e Characters? 
— A string of size 1 


e Matlab has more types than Octave, e.g. user-defined OO-classes 
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Creating a Matrix 


Simply type: 
OODaVvesl» A = [o, wy LF Dy Lacs Y, O64 =o] 


A = 
8 2 1 
>. il 4 
J 6 -5 


More on matrices, further down this tutorial. 
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Creating a Character String 


e Simply type 


octave:4> str = 'Hello World' 


Opposed to Matlab, Octave can also deal with double quotes. 
For compatibility reasons: always use single quotes 


Creating a Structure 
e Type for instance 


octave:5> data.id = 3; 
octave:6> data.timestamp = 1265.5983; 


octave: /> data.name = "sensor 1 front: 


Variables and Data Types 


Creating a Vector of Structures 


e Oh, a new measurement has arrived. Extend struct by: 


octave:8> data(2).id = 4; 
octave:9> data(2).timestamp = 1268.9613; 


octave..> data(2).name = "sensor 1 back'; 


Octave will respond with: 


data = 
{ 


1x2 struct array containing the fields: 


id 
timestamp 


name 
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Display Variables 
e Simply type its name 
octave:1> a 


a = 4 


Suppress Output 
e Addasemicolon 


octave:2> a; 


octave:3> sin(phi); 


e Applies also to function calls 
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e Variables have no permanent type. Octave/Matlab are weakly typed 
languages 


s = 3 followed by s = 'octave' isfine 


* Use who (orthe more detailed whos ) to list the currently defined 
variables. Example output: 


Variables in the current scope: 


Attr Name Size Bytes Class 
A 3x3 72 double 
a 1x1 8 double 
ans 21x1 168 double 
S 15 3 Char 


V 1x21 24 double 
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Numerical Precision 


Variables are stored as double precision numbers in IEEE floating point format 


e realmin Smallest positive floating point number: 2.23e-308 
e realmax Largest positive floating point number: 1.80e+308 
e eps Relative precision: 2.22e-16 


e These keywords are reserved and can be used in your code 
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Control Display of Float Variables 


* format short Fixed point format with 5 digits 

e format long Fixed point format with 15 digits 

e format short e Floating point format, 5 digits 

* format long e Floating point format, 15 digits 

e format short g Best of fixed or floating point with 5 digits 


(good choice) 


e format long g Best of fixed or floating point with 15 digits 


See help format for more information 
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Talking about Float Variables... 


* ceil (x) Round to smallest integer 
not less than x 


e floor (x) Round to largest integer 

not greater than x 
e round (x) Round towards nearest integer 
e fix (x) Round towards zero 


If x is a matrix, the functions are applied to each element of x. 
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Creating a Matrix 


e Simply type: 
OCtavesl> A = Lo, £r Lf Dy ~l 4 7, Op =] 


e To delimit columns, use comma or space 


e To delimit rows, use semicolon 
The following expressions are equivalent: 


A= De 2 173 =L 7 6 =p] 
A = [Ope pli op lra Oyo. 


21 


RG 


e Human-Oriented Robotics TA 1 
Matrices Prof. Kai Arras | i | "^ 


Social Robotics Lab z 
2 


FREIBU 


Creating a Matrix 


e Octave will respond with a matrix in pretty-print: 


A = 
8 2 1 
J =l 4 
+ SU 


e Alternative Example: 
octave:2» phi = pi/3; 
octave:3> R = [cos (phi) -sin(phi); sin(phi) cos(phi) |] 
| 二 
0250000 =0.86603 
0.86603 0.50000 
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Creating a Matrix from Matrices 


octave:1» A = [1 11; 22 2]; B= [33; 33]; 


e Column-wise 


octave:2> C = [A B] 
C — 
1 1 1 33 
2 2 2 3 
e Row-wise: 


octave:3> D = [A; [44 44 A4]] 
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Indexing 


Always "row before column"! 


aij = A(1i,j) Get an element 
jo XU. =) Get a row 

Se = s.) Get a column 
B= A(i:k; JIL) Get a submatrix 


Useful indexing command end : 


octave:1» data = [4 -1 35 9 11 -2]; 
octave:2> v = data(3:end) 
NV = 


os 2 Lk sm 
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The two meaning of colon : 


Wildcard to select entire matrix row or column 


A(3,:), E272) 


Defines a range in expressions like 

indices = 1:5 Returns row vector 1,2,3,4,5 

steps = 1:3:61 Returns row vector 1,4,7,...,61 

t = 0:0.01:1 Returne vector UU01L,0.02.-& 2-1 
FOX x 


start increment stop 


Useful command to define ranges: 1inspace 
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Matrices 


Assigning a Row/Column 
e All referenced elements are set to the scalar value. 


[Ll NS o Oy 
-3; 


octave:1> A = Ze ee de SS oS oS a] 


octave:2> A(3,:) = 


Adding a Row/Column 
e |f the referenced row/column doesn't exist, it's added. 


octave:3> A(4,:) = 4 
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Deleting a Row/Column 


e Assigning an empty matrix [] deletes the referenced rows or columns. 
Examples: 


octave:4> A(2,:) = [] 
A = 


4 4 4 4 4 


octave:4> A(1:2:5,:) = [] 
AS 

2 4 

2 2 

=3 =3 
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Matrices 


Get Size 


nr = size(A,1) 

nc = size(A,2) 
[nr nc] = size (A) 
l = length (A) 
numel (A) 


isempty (A) 


Octave only: 
e nr = rows (A) 
e nc = columns (A) 
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Get number of rows of A 

Get number of columns of A 
Get both (remember order) 
Get whatever is bigger 

Get number of elements in A 


Check if A is empty matrix [] 


Get number of rows of A 


Get number of columns of A 


jab 感 
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Matrices 


Matrix Operations 


B = 3*A 

C= A*B + X = 
B = A' 

B = inv (A) 

Ss = sey 

d = det (A) 

[v lambda] = 


[U S V] = svd (A) 


many many more... 


D 


e190 (A) 
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Multiply by scalar 
Add and multiply 
Transpose A 
Invert A 


Mix vectors and matrices 


Determinant of A 
Eigenvalue decomposition 


Singular value decomposition 


al 
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Vector Operations 


With x being a column vector 


es = x'*x Inner product, result is a scalar 
e X = x*x' Outer product, result is a matrix 
ee = xX*x Gives an error 


Element-Wise Operations 


es = X.+x Element-wise addition 

e p = X.*x Element-wise multiplication 
e q = x./x Element-wise division 

e e = x.^3 Element-wise power operator 
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Matrices 


Useful Vector Functions 


sum (v) 


cumsum (v) 


proc Ww.) 


cumprod (v) 


diff(v) 


mean (v) 


std(v) 
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Compute sum of elements of v 


Compute cumulative sums of 
elements of v (returns a vector) 


Compute product of elements of v 


Compute cumulative products of 
elements of v (returns a vector) 


Compute difference of subsequent 
elements [v(2)-v(1) v(3)-v(2) ...] 


Mean value of elements in v 


Standard deviation of elements 


al 
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Useful Vector Functions 


e min (v) Return smallest element in v 

* max (v) Return largest element in v 

e sort(v,'ascend') Sort in ascending order 

e sort (v, 'descend') Sort in descending order 

e find(v) Find indices of non-zero elements. 
Great in combination with vectorization 
Example: 


ivec = find(datavec == 5) 
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Special Matrices 


e A = zeros (m,n) 

e B = ones(m,n) 

e I = eye (n) 

e D = diag([a b cl) 
Just for fun 

e M = magic (n) 
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Zero matrix of size m xn 
(Often used for preallocation) 


Matrix of size m x n with all 1's 
Identity matrix of size n 


Diagonal matrix of size 3 x 3 
with a,b,c in the main diagonal 


Magic square matrix of size n x n. 
(All rows, columns sum up to same number) 
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Matrices 


Random Matrices and Vectors 


e R 


rand (m,n) 


randn (m,n) 


randperm (ín) 
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Matrix with m x n uniformly 
distributed random numbers 
from interval [O..1] 


Matrix with m x n normally 
distributed random numbers 
with zero mean, unit variance 


Row vector with a random permutation 
of the numbers 1 ton 
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Multi-Dimensional Matrices 

Matrices can have more than two dimensions. 

e Create a 3-dimensional matrix by typing, e.g., 
octave:1> A = ones(2,5,2) 


Octave will respond by 


1 1 1 1 1 

il 1 1 il 1 
ans(:,:,2) = 

i 1 1 il 1 

1 1 1 1 1 


35 


e Human-Oriented Robotics ia 
Matrices Prof. Kai Arras | i | T 


Social Robotics Lab 


FREIBURG 


UNI 


Multi-Dimensional Matrices 


e All operations to create, index, add, assign, delete and get size 
apply in the same fashion 


Examples: 
e [nn 1] = s1ze(A) 
e A = ones (m,n, 1) 


e m = min(min(min(A) )) 


* aijk = A(i,j,k) 
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Matrix Massage 


Matrix operations that have no mathematical meaning. Useful for 
manipulating data which is organized in matrices 


e reshape (A,m, n) Change size of matrix A to have dimension 
m x n. An error results of A does not have 
m x n elements 


e circshift(A,[m n]) Shift elements of A m times in row 
dimension and m times in column 
dimension. Has no mathematical meaning 


e shiftdim(A,n) Shift the dimension of A by n. 
Generalizes transpose for multi- 
dimensional matrices 
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Matrix Massage 


。 fliplr (A) 


e flipud (A) 


e flipdim (A, dim) 


e rot90 (A) 
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Reverses the order of columns of matrix A 
in left/right-direction. Rows are not changed 


Reverses the order of rows of matrix A 
in up/down-direction. Columns are not 
changed 


Flip matrix A along dimension dim. Typically 
for multi-dimensional matrices 


90 degree counterclockwise rotation of 
matrix A. This is not the transpose of A 
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Matrix Massage Example 


Let P = [x1; y1; x2; y2; ...] be a 2nx1 column vector of n (x,y)-pairs. 
Make it a column vector of (x,y,theta)-tuples with all theta being pi/2 


e Make P ita 2 xn matrix 


octave:1» P = reshape(P,2,numel (P)/2) ; 


e Add a third row, assign pi/2 
octave:2> P(3,:) = pi/2; 


e Reshape it to be a 3n x 1 column vector 


octave:3> P = reshape(P,numel(P),1); 
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Most Often Used Commands 


strcat Concatenate strings 

int2str Convert integer to a string 

num2str Convert floating point numbers to a string 
Sorin. Write formatted data to a string. 


Same as C/C++ fprintf for strings 


Example 

S = strcat('At step ',int2str(k),', p = ',num2str(p,4)) 
Given that strings are matrices of characters, this is equivalent to 
S = ['At step ' int2str(k) ', p = ' num2str (p, 4) ] 
Octave responds with 


S = At step 56, p = 0.142 
40 
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Octave/Matlab has virtually all common string and parsing functions 


e You are encouraged to browse through the list of commands or simply 
type help command: 


strcmp, strncmp, strmatch, char, ischar, findstr, 
strrind, sStrZdouble, strznum, numZstr, Strvcat, 
strtrim, strtok, upper, lower,... 


and many more... 
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Plotting in 2D 
e plot (x, cos (x)) Display x,y-plot 


Creates automatically a figure window. Octave uses gnuplot to handle 
graphics. 


e figure (n) Create figure window 'n' 


If the figure window already exists, brings it into the foreground 
(= makes it the current figure) 


e figure Create new figure window with 
identifier incremented by 1 
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Several Plots 
e Series of x,y-pairs: plot (xl,yl,x2,y2,...) 


e.g. plot(x,cos(x),x,sin(x),Xx,x.^2) 


e Add legend to plot: command legend 


legend('cos(x),'sin(x)','x^2') 


e Alternatively, hold on does the same job: 


octave:l» hold on; plot(x,cos(x)); 
octave:2> plot(x,sin(x)); 


Octdveroc pDhobtxx. 2)? 
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Plotting 


Frequent Commands 


cli 
hold on 


grid on 


grid off 


title ('My Plot') 
xlabel('time') 


ylabel ('prob') 
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Clear figure 


Hold axes. Don't replace plot with 
new plot, superimpose plots 


Add grid lines 


Remove grid lines 


Set title of figure window 
Set label of x-axis 


Set label of y-axis 


IIA- 
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Controlling Axes 


e axis equal Set equal scales for x-/y-axes 
(Use it!) 

* axis square Force a square aspect ratio 

* axis tight Set axes to the limits of the data 

* à — axis Return current axis limits 


[xmin xmax ymin ymax] 
e axis([-1 1 2.5 5]) Set axis limits (freeze axes) 


e axis off Turn off tic marks 


e box on Adds a box to the current axes 


e box off Removes box 
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Controlling Plot Styles 


e In plot (x, cos (x), 'r+') the format expression 'r+' means 
red cross 


e There are a number of line styles and colors, 
see help plot 


Example: 
octave:l> x = linspace(0,2*p1,100) ; 


Sotave > slot (x cos (x), TEE, x sinz); bat) 


produces this plot: 
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PLOCI; COS) r Er usc usu o6] v OX]? 
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e Adjusting the axes 


OCLave:o» axis ( [0 2*pi -1 1]) 


(try also axis tight) 


e Adding a legend, labels and a title 


octave:4» 
Legend('eostx)','ssnix)'," Locatron'.,"'Southwest') 


octave:5> title('Trigonometric Functions') 
octave:6» xlabel('x') 


octave: 7> ylabel('y') 
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Uhm..., don't like it. Let’s start over... 


octave:l> clf; 


e Controlling Color and Marker Size 


OCtadve.2> lot (x, COS(X) (EE p< 7S nr X 7 ces 
'Color',[1 .4 .8], 'MarkerSize',2) 


octave:3> axis tight 


e Adding Text 
octave:4> text(1,-0.5,'cos(\phi)') 
ctave5> ext, 0.5, exa Com) EN 
Note the LaTeX syntax! 
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sin(@) £f 
na 
na 
十 
十 
十 
十 
十 
十 
十 
+ 
+ 
十 
十 
十 
Bin 
E 
-1 | | | eene | 
0 1 3 4 5 6 


DiLOL (x, COS (x) y Prt" sr), FX, Color"; [1 .4 0l 'Marker ire", 2) 
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Plotting 


Yepp, | like it... Get hardcopy! 


Exporting Figures 

e print -deps myPicBW.eps 

e print -depsc myPic.eps 

e print -djpeg -r80 myPic.Jjpg 
e print -dpng -r100 myPic.png 
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Export B/W .eps file 
Export color .eps file 
Export .jpg in 80 ppi 
Export .png in 100 ppi 


See help print for more devices including specialized ones for Latex 


。 print can also be called as a function. 


Then it takes arguments and options as a comma-separated list. 


DroneTt'e-—ng' Secus Spo mts 
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This tutorial cannot cover the large variety of graphics commands in 
Octave/Matlab 


e You are encouraged to browse through the list of commands or simply 
type help command: 


hist, bar, pie, area, fill, contour, quiver, 
Scatter, compass, rose, semilogx, loglog, stem, 
stairs, image, imagesc 


and many more! 
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Plotting in 3D 

e plot3 Plot lines and points in 3d 
e mesh 3D mesh surface plot 

e surf 3D colored surface plot 


Most 2d plot commands have a 3D sibling. Check out, for example, 


bar3, pie3, fill3, contour3, quiver3, 
Scatter3, stem3 


Let us look at some examples... 
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Example: plot 


Morse Signal Analysis 


$ Load data 300 


load MDdata xdata distl dist2 dist3 O Stress 
+ Sammon Mapping 


250 A Squared Stress 
% Plot the first set of data in blue 


figure; hold on; 

plot(xdata, distl, 'bo'); 
plot(xdata, dist2, ‘r+'); 
plot(xdata, dist3, 'g^'); 


200 


Distances 
em 
e 


$ Add title, axis labels, legend 100 
title('Morse Signal Analysis'); 
xlabel('Dissimilarities'); 50 


ylabel('Distances'); 


legend(('Stress', 'Sammon Mapping', 0 (Č 
0 20 40 60 OÙ 100 120 140 160 180 200 


‘Squared Stress'),'Location','NorthWest'); Dissimilarities 
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Example: plot3 


% Load data 
load SpectraData massc time spectra; 


Extracted Spectra Subset 


$ Create the 3D plot 

figure; 

plot3(massc, time, spectra); 
box on; 


$ Set viewing angle and axis limits 
view(26, 42); 
axis([500 900 0 22 0 4e8]); 


lon Spectra 


$ Add title and axis labels 
xlabel('Mass/Charge (M/Z)'); 
ylabel('Time'); 

zlabel('Ion Spectra'); 
title('Extracted Spectra Subset'); 


Mass/Charge (M/Z) 
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Example: ezplot 


Lemniscate Function 
G2 + 22 - x2 + ÿ2 


$ Create the plot 

figure; 

ezplot('(x^2 + y^2)^2 = x^2 + y^2' ,... 
[-1.1, 1.1], [-1.1, 1.1]); 


% Add a multi-line title 
title({'Lemniscate Function';... 


"(u(*-" 2 b v2) 2 «2 F y 2 |}; 


Note: the special character ... attheendofaline 
continues the current function on the next line 
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Example: bar 


$ Load data 4 


x 10 Childhood diseases by month 
load Datafile measles mumps chickenpox; 
BE Measles 
o E Mumps 
% Create a stacked bar chart bar EB Chicken "1 


figure; 
bar(1:12, [measles mumps chickenpox],... 
0.5, 'stack'); 


$ Adjust the axis limits 
axis([0 13 0 100000]); 


Cases (in thousands) 


Q 


% Add title, axis labels, legend 
title('Childhood diseases by month'); 
xlabel('Month'); 

ylabel('Cases (in thousands)'); 


legend('Measles', 'Mumps', 'Chicken pox'); 


Month 
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Example: bar3 


$ Load monthly temperature data 


O 


Boston Monthly Temperatures 1900-2000 


load MonthlyTemps temperatures months years; 


$ Create the 3D bar chart 
figure; 
bar3(temperatures); 


axis([0 13 0 12 0 80]); 


$ Add title and axis labels 
title('Boston Monthly Temps 1900-2000'); 
xlabel('Month'); 

ylabel('Year'); 

zlabel('Temperature'); 


Temperature 


$ Change the x and y axis tick labels 
set(gca, 'XTickLabel', months); 
set(gca, 'YTickLabel', years); 
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Example: polar 


abs(sin(2t)*cos(2t)) 
90 05 
$ Create data for the function - 
t = 0:0.01:2*pi; 
r - abs(sin(2*t).*cos(2*t)); 


oo 


Create a polar plot using polar 
figure; 


polar(t, r); 


$ Add a title 


title('abs(sin(2t)*cos(2t))'); 270 
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Example: scatter3 


Ozone Levels 


% Load data 14 ppm 


o 


load OzoneData ozoneidx temp wind rad; 
12 ppm 
$ Create a 3D scatter plot 
figure; 10 ppm 
scatter3(temp, wind, rad, 30, ... 
ozoneidx, 'filled'); 


view(-34, 14); 


8 ppm 


Solar Radiation 


6 ppm 
$ Add title and axis labels 
title('Ozone Levels'); 

xlabel( 'Temperature' ); "pom 
ylabel('Wind Speed'); 


zlabel('Radiation'); 


2 ppm 


Wind Speed Temperature 


$ Add a colorbar with tick labels 
colorbar('location', 'EastOutside', 'YTickLabel',... 
{'2 ppm', ‘4 ppm, ‘6 ppm', '8 ppm, '10 ppm', '12 ppm', '14 ppm')); 


For individually colored points, use scatter instead of plot in a for-loop! 
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Example: surfc 


oo 


Create a grid of x and y data 
-10:0.5:10; Normal ion 
-10:0.5:10; ae 
X, Y] = meshgrid(x, y); 


y 
x 
[ 


$ Create the function Z = f(X,Y) 
Z = sin(sqrt(X.^2*4Y.^2))./sqrt(X.^2*Y.^2); 


$ Create a surface contour plot 
figure; 

surfc(X, Y, 2); 

view(-38, 18); 


$ Add title and axis labels 
title('Normal Response'); 
xlabel('x'); 

ylabel('y'); 

zlabel('z'); 
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e Start, quit, getting help 
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e Octave and Matlab in practice Matlab 


e librobotics 
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Programming in Octave/Matlab is super easy 


e But keep in mind: indexing is one-based, i.e. 
Indices start with 1 !!! 


OcqcdVveslo vy = 1210 
octave:2> v(0) 


error: subscript indices must be either positive integers or 
logicals 


e Octave/Matlab is case-sensitive 


Text Editors 
e Use an editor with m-file syntax highlighting/coloring 
e Matlab has its own IDE 
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Control Structures 
* if Statement 


it GOUnOdIUtiofn, 
then-body; 
elseif condition, 
elseif-body; 
else 
else-body; 


end 


e The else and elseif clauses are optional 


e Any number of elseif clauses may exist 
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Control Structures 
* switch Statement 


Switch expression 
case label 

command-list; 
case label 


command-list; 
otherwise 


command-list; 


end 


e Any number of case labels are allowed 
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Control Structures 


while Statement 


while condition, 
body; 


end 


for statement 


for var = expression, 
Dody; 


end 
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Interrupting and Continuing Loops 


e break 


Jumps out of the innermost for or while loop that encloses it 


e continue 


Used only inside for or while loops. It skips over the rest of the loop 
body, causing the next cycle to begin. Use with care 
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Increment Operators (Octave only!) 


Increment operators increase or decrease the value of a variable by 1 


T4 


Increment scalar i by 1 
Decrement scalar i by 1 
Increment all elements of matrix A by 1 


Decrement all elements of vector v by 1 


There are the C/C++ equivalent operators ++i, --A 
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Comparison Operators 


e All of comparison operators return a logical value of 1 if the comparison 
is true or a logical value of 0 if it is false 


1 == 6, condl = (d > theta) 


e For the matrix-to-matrix case, the comparison is made on an 
element-by-element basis 


[1 2; 3 4] == [1 3; 2 4] returns [1 Os O 1] 


e For the matrix-to-scalar case, the scalar is compared to each element 
in turn 


[1 2; 3 4] == 2 returns [O 1; O QO] 
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Comparison Operators 


* any(v) Returns 1 if any element of 
vector v is non-zero (e.g. 1) 


e all(v) Returns 1 if all elements in 
vector v are non-zero (e.g. 1) 


For matrices, any and all return a row vector with elements corresponding 
to the columns of the matrix 


* any (any (C)) Returns 1 if any element of 
matrix C is non-zero (e.g. 1) 


e all(all(C)) Returns 1 if all elements in 
matrix C are non-zero (e.g. 1) 
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Relational Operators 


e X< y 

s x y 
ex == y 
" XxX 2 y 
e X> y 

e x ~= y 
e x l= y 
e X <> y 
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True if x is less than y 


True if x is less than or equal to y 


True if x is equal to y 
True if x is greater than or equal to y 
True if x is greater than y 


True if x is not equal to y 


True if x is not equal to y (Octave only) 


True if x is not equal to y (Octave only) 


al 
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Boolean Expressions 


e Bl & B2 Element-wise logical and 

e B1 | B2 Element-wise logical or 

e ~B Element-wise logical not 

e !B Element-wise logical not (Octave only) 


Short-circuit operations: evaluate expression only as long as needed 
(more efficient) 


e Bl && B2 Short-circuit logical and 


e B1 || B2 Short-circuit logical or 
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Recommended Naming Conventions 
e Functions: underscore-separated or lowercase notation 


Examples: drawrobot.m, calcprobability.m, 
Intersect Ae OBL. 


e Scripts: UpperCamelCase 


Examples: LocalizeRobot.m, MatchScan.m 


e Matlab/Octave commands are all in lowercase notation 
(no underscores, no dashes) 


Examples: continue,int2str, isnumeric 
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Functions 


Octave/Matlab programs can often be simplified and structured by defining 
functions. Functions are typically defined in external files, and can be called 
just like built-in functions 


e Inits simplest form, the definition of a function looks like this: 


Function name 
body 


end 


e It is recommended to define one function per file 


e These files are called m-file or .m-file 
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Passing Parameters to/from Functions 


Simply write 
function [ret-var] = name(arg-list) 
body 


end 


arg-list is acomma-separated list of 
input arguments arg1, arg2, ..., argn 


ret-var is a comma-separated list of 
output arguments. Note that ret-var is a vector enclosed in square 
brackets [arg1, arg2, ..., argm]. 
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Functions and Scripts 


Examples Please: 


function [mu sigma] = 
mu = mean (data); 
Sigma = std(data); 

end 

function [haspeaks i] 
indices = find(data 


1f isempty (indices), 


haspeaks = 0; 1 = 
else 

haspeaks = 1; i = 
end 


end 
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calcmoments (data) 


findfirstpeak (data, thresh) 


> thresh); 


Oe 


indices (1); 
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Local Variables, Variable Number of Arguments 


e Of course, all variables defined within the body of the function are 
local variables 


e varargin Collects all input argument in a cell 
array. Get them with vararginti} 


* varargout Collects all output argument in a cell 
array. Get them with varargoutii} 


e nargin Get the number of input args 


e nargout Get the number of output args 


e See help varargin, help varargout for details 
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Functions and their m-File 


e When putting a function into an m-file, the name of that file must be 
the same than the function name plus the .m extension 
Examples: calcmoments.m, findfirstpeak.m 


e To call a function, type its name without the .m extension. 
Example: 


[bool 1] = findfirstpeak(myreadings, 0.3); 


e Comments in Octave/Matlab start with % . Use them a lot! 
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Scripts 


The second type of m-files is called script. Again, Octave/Matlab scripts 
are text files with an .m extension 


Scripts contain executable code. They are basically the "main" programs 


Execute a script by typing its name without the .m extension 
Example: octave:1> LocalizeRobot 


Again, comments in Octave/Matlab start with 5. 
(I can't repeat this often enough ;-) 
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Document your Function/Script 


e You can add a help text to your own functions or scripts that then 
appears on help command 


e The first block of comment lines in the beginning of an m-file is defined to 
be help text. Example: 


SNORMANGLE Put angle into a two-pi interval. 


[9] 


$ AN = NORMANGLE(A,MIN) puts angle A into the interval 


% MIN MINE “pals If A avs int, Ink 1S returned. 


% v.1.0, Dec. 2003, Kai Arras. 


function an = normangle(a,mina) ; 
if a < Inf, 


issu] 


help text 
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Setting Paths 


path 


addpath('dir') 


rmpath('dir') 


savepath 
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Print search path list 


Prepend the specified directory 
to the path list 


Remove the specified directory 
from the path list 


Save the current path list 
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Save Variables 


After a complex or lengthy computation, it is recommended to 
save variables on the disk 


* Save my varssmat 
Saves all current variables into file my vars.mat 


e save results.mat resultdata X Y 
Saves variables resultdata, X and Y in file results.mat 


e save ... -ascii 
Saves variables in ASCII format 


e Seve do» Mat 
Saves variables in binary MAT format 
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Load Variables 


The corresponding command is 1oad 


e load my vars.mat 


Retrieves all variables from the file my. vars.mat 


e load results.mat X Y 


Retrieves only X and Y from the file results.mat 


An ASCII file that contains numbers in a row/column format (columns 
separated by spaces or commas, rows separated by new lines) can be 
simply read in by 


e A = load('data.txt') 


Matrix A will then contain the data 
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Open, Write, Close Files 


e fopen Open or create file for writing/reading 
e fclose Close file 
& Lorin i Write formatted data to file. C/C++ 


format syntax 


Example: 
v = randn (1000,1); 
fid = fopen('gauss.txt','w'); 
for a = l:length(v)y 
fprintf (fid; CJ Jp 7371) yy 
end 
fclose(fid); 
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Attention, Popular Bug 


e If your program writes to and reads from files, floating point 
precision of fprintf is crucial! 


e Be sure to always write floating point numbers into files using the 
appropriate precision 


e [n the above example, with format definition '%7.4£\n', this file 
will be a very poor source of Gaussian random numbers 
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Reading Files (more advanced stuff) 


e textread Read formatted data from text file 
e fscanf Read formatted data from text file 
e fgetl Read line from file 

e fread Read binary data file 


Read/write images 
e imread Read image from file (many formats) 


e imwrite Write image to file (many formats) 
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Cleaning Up 
* clear A Clear variable A 
* clear frame* Clear all variables whose names 
start with frame, e.g. frame001, frames 
* clear Clear all variables 
* clear all Clear everything: variables, 
globals, functions, links, etc. 
e close Close foreground figure window 
* close all Close all open figure windows 


e cle Clear command window (shell) 
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Displaying (Pretty) Messages 
e disp (A) Display matrix A without printing 
the matrix name 


e disp(str) Display string str without printing 
the string name 


Example: when typing 
octave:l> disp('done') 
Octave will print 


done 
instead of 
ans = done 


from sprintf ('done') or 'done' 
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Command History 


e Navigate up and down the command history using the up/down 
arrow keys 


e The command history is start-letter sensitive. Type one or more letters 
and use the arrow keys to navigate up and down the history of 
commands that start with the letters you typed 


Tab completion 


。 Octave/Matlab have tab completion. Type some letters followed by tab 
to get a list of all commands that start with the letters you typed 
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Built-in Unix Commands 


* pwd Display current working directory 
e ls List directory. See also dir 

e cd Change directory 

e mkdir Make new directory 

e rmdir Delete directory 


Related Commands 


e movefile Move file 


e copyfile Copy file 
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Random Seeds 


* rand and randn obtain their initial seeds from the system clock 


e To generate repeatable sequences of random numbers, set the random 
generator seeds manually 


To set the random seeds: 
e rand('seed',val) Set seed to scalar integer value val 


e randn('seed',val) Set seed to scalar integer value val 
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Useful Stuff in Practice 
We will cover: 


1. Generating output from a C/C++/Python/Java/... program in Matlab 
syntax, e.g. using Octave/Matlab as a visualizer front-end 


2. Making animations (without Matlab’s movie function) 


3. Calling unix/dos functions from within Octave/Matlab 
programs 


4. Increasing speed through vectorization and preallocation 
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Writing Files in Matlab Syntax 


Octave/Matlab are very powerful visualization tools 


Regular languages such as C/C++/Python/Java/etc. have some support 
for graphical output but in comparison their libraries are not as flexible, 
powerful and easy-to-use than Octave/Matlab 


So, how can we combine the advantages? 


For testing or developing an algorithm in C/C++/Python/Java/etc., it is 
typically necessary to plot many variables, visualize intermediate and 
final results or make animations. Instead of writing complex 
visualizations in those languages, use Octave/Matlab as visualizer front- 
end 


Drawback: not real-time (can be made quasi real-time) 
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Writing Files in Matlab Syntax 


e Data written into plain text file in matrix format. 
Example: 


LIJIDOrec CAL XL 


Oe 22756 D D S W295 OS 0.647680 DEUS D 
Qe 0679 0-414615 0 70472 DES ADEM Qu Sk om 
Ore Cols OO 0062014 Dos uod GES Zils SZ OFS 70605 
0.940364 ROSES GES d oL Ve 2306906 OS 604 7 20 
ORs Sic 0-872248 Dados 0.042745 02228580 


e Read in using the command load. 
Example: A = load('filtered readings.txt'); 
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PlotFilteredReadings.m 


e File may also contain Matlab code snippets. Example: 


A = 


lG 


trigue ly Ce c Bold on; 
lotl: sIizelA AE 


[ 


SRS SRS 


-1792238 
2290908 
CUS UE 
- 940364 


SES => 


02-029 
S41 A615 
20/65 914 
319/515 5E S 


© O ere 


-997623 
210472 
es od 
6 Suet Ik 


Oo 0 © © 


-647080 
SON D 
25 7155.57. 
20106 S016 


e Must have the .m extension. It's a script. 


O O © O 


-198282 
-043092 
Fie ols 
.804710 


e Simply execute by typing PlotFilteredReadings 
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Making Animations 

e Matlab has commands such as getframe and movie to make 
animations from plots 

e Octave, being free of charge, does not (yet) support these commands 


e Never mind! Here is a pretty obvious way to make movies: 


— Export plots to a directory (e.g. “frames”) using print 
from within a loop. Then compose frames to a movie using 
tools such as ImageMagick or Quicktime Pro. 
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Making Animations. Example: 


e Let data.txt contain data in matrix format, we want to plot each 
column and save it as a frame 


A = load('data.txt'); 

[nrow ncol] = size(A); 

figure (1); 

for 1 = 1:ncol， 
DLOEt (LL nNrOw; ACs me 
fname = sprintf('frames/frame$04d.png',i); 
Prine 70009; UUSee) 


end 


e Problem: axis limits change for each plot/frame. 
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Making Animations. Example: 


e To freeze the axes over the entire animation, use the command 
axis([xmin xmax ymin ymax]) after the plot command 


A = load('data.txt'); 
[nrow ncol] = size(A); 
figure (1); 

for 1 = 1:ncol， 


Clee (Lent Ow, 202 2) 


axis([l1 nrow min(min(A)) max(max(A))]); 
fname = sprintf('frames/frame£04d.png',i); 


printi'—dpug"',"-:100'*',47nsme); 
end 
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Calling unix/dos Functions 


e ForUnix/Linux/Mac OS X systems, there is the command unix to 
execute system commands and return the result. 
Examples: 
unix('ls -al') 
unc c EES ~ TE DELL 
Unix (" ./myprogram’ ) 


e For Windows PCs, there is the equivalent command dos. 
e These commands allow for powerful and handy combinations with 
other programs or system commands 


e Can help to accelerate edit-compile-run cycles or edit-compile-run- 
visualize cycles in particular when Octave/Matlab is used as a visualizer 
front-end 
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Speed! 


e The low execution speed of Octave/Matlab programs is commonly 
recognized to be their most important shortcoming 


。 Mostly your program is slow, not the built-in functions! 


e This brings us to the following guidelines 


e For-loops are evil 


e Vectorization is good 


Preallocation is good 


Prefer struct of arrays over arrays of struct 


e Advanced topics (not covered here): Matlab compiler, linking C/C++, 
Fortran code from Matlab programs (mex files), parallel computing, etc. o 
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Speed: Vectorization 
e Given phi = linspace(0,2*pi,100000); 


The code 
for 1 = 1:length (phi), 
sinphi(i) = sin(phi(i)); 
end; 
is significantly slower than simply 


Sinphi = sin(phi); 


e All built-in commands are vectorized, i.e. allow vector arguments 


e You have to (and will) learn to think vectorized! 
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Speed: Preallocation 


e Ifa for- or while-loop cannot be avoided, do not grow data structures in 
the loop, preallocate them if you can. 
Instead of, for example 


for i = 1:100, 
A(1,:) = rand(1,50); 
end; 


write 


A = zeros(100,50); $ preallocate matrix 
for i = 1:100, 

A(i,:) = rand(1,50); 
end; 
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Speed: Structure of Arrays 
e Always prefer a struct of arrays over a array of structs (called plane 
organization vs. element-by-element organization) 


e It requires significantly less memory and has a corresponding 
speed benefit 


e Structure of arrays 
data.x = linspace(0,2*p1,100); 
data.y = sin(data.x); 

e Array of structure 


people(l).name = 'Polly J Harvey; 
people(l).age = 29; 


people (1000) .name = "Big Lebowski'; 
people (1000) .age = 35; 
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e Overview 


e Start, quit, getting help 


e Variables and data types 
e Matrices GNU Octave 
e Plotting 

e Programming 

e Functions and scripts 
e Files I/O 

e Misc 


e Octave and Matlab in practice Matlab 


e librobotics 
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e librobotics is a small library with frequently used Octave/Matlab 
functions in Robotics, especially for visualization 


chi2invtable.m drawrawdata.m JZ2COmp m 
compound.m drawreference.m jinv.m 
diffangle.m drawrobot.m mahalanobis.m 
drawarrow.m drawroundedrect.m meanwm.m 
drawellipse.m drawtransform.m normangle.m 
drawlabel.m 1 COMPOUNCL.TM 

drawprobellipse.m 1 VCOMp i 


e Download from SRL Homepage: 
srl.informatik.uni-freiburg.de/downloads 
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e Command drawreference.m 
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e Command drawroundedrect.m 
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e Command drawarrow.m 
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« Command drawprobellipse.m 


115 


o o Human-Oriented Paois or 
| i b robotics Prof. Kai Arras | i | 


Social Robotics Lab 


e Command drawtransform.m 
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e Command drawrobot.m 


117 


| e b b e Human-Oriented Robotics ie 
Prof. Kai A 3 
Ibrobotics es (ll 1d 
Da 


e Example Figure 
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e All commands are fully documented, 
justtype help command. 


e The command chi2invtable.m returns values of the cumulative 
chi-square distribution, often used for gating and hypothesis testing. 
It replaces the chi2inv function from the Matlab statistics toolbox 
- a costly addition to Matlab - and is also much faster 


e librobotics is compatible with both, Matlab and Octave 


e |t's open source, feel free to distribute and extend 


e Link: http://srl.informatik.uni-freiburg.de/downloads 
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Full Octave online documentation 


e http://www.octave.org 
> Support 
> Documentation 
> Reference manual in HTML or pdf (800 pages) 


e Directly: www.gnu.org/software/octave/octave.pdf (Oct 2013) 


Full Matlab online documentation: 


e http://www.mathworks.com 
> Products & Services 
> MATLAB 
> Documentation 


e Directly: http://www.mathworks.com/help/matlab/index.html (Oct 2013) 
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Thanks and Enjoy! 


Kai Arras, Social Robotics Lab 
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